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Abstract 

Basing upon dual analytic models, we present arguments in favor of 
the parametrization of generalized parton distributions (GPD) in the form ~ 
{%/where a(t) = a(t) — a(0) is the nonlinear part of the Regge trajectory 
and go is a parameter, go > 1. For linear trajectories it reduces to earlier proposals. 
We compare the calculated moments of these GPD with the experimental data on 
form factors and find that the effects from the nonlinearity of Regge trajectories are 
large. By Fourier-transforming the obtained GPD, we access the spatial distribution of 
protons in the transverse plane. The relation between dual amplitudes with Mandel¬ 
stam analyticity and composite models in the infinite momentum frame is discussed, 
the integration variable in dual models being associated with the quark longitudinal 
momentum fraction x in the nucleon. 
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1 Introduction 


Generalized parton distributions (GPD) [TJ -2, fcf combine our knowledge about the one¬ 
dimensional parton distribution in the longitudinal momentum with the impact-parameter, 
or transverse distribution of matter in a hadron or nucleus. It is an ambitious program 
to access the spatial distribution of partons in the transverse plane and thus to provide a 
3-dimensional picture of the nucleon (nucleus) 7|8]. This program involves various 

approaches, including perturbative QCD, Regge poles, lattice calculations etc. (see Ref. [Hj 
for reviews). The main problem is that, while the partonic subprocess can be calculated 
perturbatively, the calculation of GPDs require non-perturbative methods. GPDs enter in 
hard exclusive processes, such as deeply virtual Compton scattering (DVCS); however, they 
cannot be measured directly but instead appear in convolution integrals, that cannot be 
easily converted. Hence the strategy is to guess the GPD, based on various theoretical 
constraints, and then compare it with the data. In the first approximation, the GPD is 
proportional to the imaginary part of a DVCS amplitude, therefore, as discussed in pm. 
the knowledge (or experimental reconstruction) of the DVCS amplitude may partly resolve 
the problem, provided the phase of the DVCS amplitude is also known. In other words, a 
GPD can be viewed as the imaginary part of an antiquark-nucleon scattering amplitude, or 
a quark-nucleon amplitude in the u channel. 

Alternatively, one can extract mm, still in a model-dependent way, the nontrivial 
interplay between the x and t dependence of GPD from light-cone wave functions 

H{x , £ = 0 ,t) = I d 2 k_L^*(^, kj_)^(a;, k ± ) + (1 - x)q ± ), 

where ^(x, kj_) is a 2-particle wave function (see, e.g., [131) and t = qj_. 

In two recent papers HUES various forms of GPD for £ = 0 were tested against the 
experimental data on the related form factors. The agreement with the data in Ref. |[3| 
is impressive; in Ref. H3 the spatial distribution of partons in the transverse plane was 
also calculated. We pursue the approach of Refs. nu m by bringing more arguments 
coming from duality in favor of the parametrization for H(x,t) used in [31 and exploring 
how analyticity and unitarity affect the t dependence of GPD, the observable form factors 
and the calculated distribution of partons in the longitudinal and transverse planes. The 
effects are very large. Similarly to papers hu m , we limit ourselves the case of vanishing 
skewedness, £ = 0. 

Regge trajectories play a key role in this analyses. Actually, there are two groups of 
trajectories in the problem: one, the p, u ; etc. trajectories exchanged in the valence quark 
distribution function H(x,t) (or, equivalently, the imaginary part of the qp scattering am¬ 
plitude). Less evident are the characteristics of the corresponding trajectory in ’’magnetic” 
densities E(x,t) for they cannot be expressed in terms of any known parton distribution. 
In Ref. [T5| the p trajectory, a p (t ) = 0.48 + 0.88t, was fitted to the masses of p( 770) and 
p 3 (1690) and a u {t) = 0.42 + 0.95t to o;(782) and cn(1690). In [3|, instead, the slopes of the 
trajectories in H(x,t ) and £(x,t) were fitted to the data on form factors and are equal to 
1.098 GeV 2 and 1.158 GeV 2 , respectively. The relevant intercepts are contained in the quark 
distributions, as discussed in Sec. 5. The resulting GPDs and related observables are very 
sensitive to the above parameters (see Fig. 1 in Ref. [T4])- Even more sensitive are they to 


any deviation from linear trajectories, as shown in Sec. 6. 

One can start either from trajectories fitted to resonances and scattering data 
(parametrizations of non-linear complex mesonic trajectories fitting the spectra of resonances 
as well scattering data exist in the literature main!) or treat them as ’’effective” ones, to 
be fitted to the data on form factors. Our strategy here is to start from trajectories close to 
those in |T2| (fitted to form factors) and then look for the effects coming from the deviations 
from linearity. 

The paper is organized as follows. In Sec. 2, following earlier publications, we show how 
dual models with Mandelstam analyticity can be related to deep inelastic scattering and 
GPD. Regge trajectories satisfying the analyticity and unitarity constrains are introduced 
in Secs. 2 and 3, where their role in the calculation of form factors is also discussed. The 
relation between GPD and form factors is discussed of Sec. 4. A particular model of GPD 
with its t dependence determined by dual models with Mandelstam analyticity, introduced in 
Sec. 2, is discussed in Sec. 5. Numerical calculations (form factors and parton distributions) 
are presented in Sec. 6. Our (temporary) conclusions and a discussion can be found in 
Section 7. 


2 Non-linear Regge trajectories 
in dual and composite models 


Dual Amplitudes with Mandelstam analyticity (DAMA) were suggested (see jTHj and earlier 
references therein) as a way to solve the manifestly non-unitarity of narrow-resonance dual 
models ns). The (u, t ) term of the crossing-symmetric DAMA is 


D(u, t ) 



— x)) — 1 


1 — X 
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—a(ux) — l 




(i) 


where u and t are the Mandelstam variables, and rp, g 2 are parameters, g \, g 2 > 1. In what 
follows we set, for simplicity, gi — g 2 — go- Similar expressions are valid for the (si) and (su) 
terms. They are not unique since the integrand of m can be multiplied by functions of the 
type f(t( 1 — x))f(ux). Furthermore, the powers in the integrand can be shifted by integers 
determined by the quantum numbers of a particular reaction and relevant exchanges. 


The functions a(y,x )), y = s, t, u, called in jTSj homotopies, map the physical Regge 
trajectories a(y ) onto linear functions a + by. Contrary to the narrow-resonance (=linear 
trajectories) Veneziano amplitude |19] . applicable only to soft collisions of extended objects 
(strings), and decreasing exponentially at any scattering angle, DAMA does not only allow 
for, but even requires the use of nonlinear Regge trajectories. It will be convenient to write 
a(y ) = a(0) + a(y), where a(y) will denote the nonlinear part of the trajectory. 

For |w| —> oo and fixed t, DAMA is Regge-behaved 

D(u,t) ~ gl(-gou) a[t) {G(t) + ...], 


G(t) = 


dyy v y a(t) \ 


where 



provided PH 


< const, 

— XX) 

which is equivalent to saying that the real part of the trajectory is bounded 1 . Compatibility 
with the wide-angle scaling behavior of the amplitude, typical of point-like constituents, 
lowers this growth to a logarithm. Examples will be presented below. 

The pole structure of DAMA 


a(u ) 
y/ulnu 


D(u,t) ~ 

1=0 


C n -l(t) 

[n — a(u)] l+1 ’ 


( 2 ) 


where C n -i(t ) is the residue, whose form is fixed by the dual amplitude (see (IHJ), is similar 
to that of the Veneziano model except that multiple poles appear on daughter levels ra- 
The pole term (2) in DAMA, comprising a whole sequence of resonances lying on a complex 
trajectory a(u), is a generalization of the Breit-Wigner formula. Such a ’’reggeized” Breit- 
Wigner model has little practical use in the case of linear trajectories, resulting in an infinite 
sequence of poles, but it becomes a powerful tool in case of complex trajectories with a 
limited real part and hence a limited number of resonances. 

The threshold behavior of DAMA satisfying the unitarity constrains [221 

D(u, t ) ~ y/t — to[const + ln(l — to/t)\, 

t—*to 


is correlated with that of the trajectories [T%] , 

9fa(t) ~ (t - f 0 ) Q(to)+1/2 - 

t—rta 


For a light threshold this is close to the square-root behaviour to be used below. 

A simple model, compatible both with to the above threshold behavior and with the 
Regge asymptotics (or polynomial boundedness) of the amplitude, yet fitting the data on 
resonances spectra, can be made of a sum of square roots 2 (the assignment of the signs is 
uniquely determined by the requirement of positivity of the imaginary part etc., see |2Tj ) 

m 


a(t) = q(0) - ^ QitVf T-t- VU)- (3) 

i 

Linear trajectories appear as the limiting case of an infinitesimally heavy threshold t\ —> oc 
in 

<a(t) = ol o — cxi (y/ti — t — y/t\) 

with fixed forward slope of the trajectory ol = The limit of the infinitely rising linear 
trajectory is associated with a hadronic string, while the finite value of t t can be interpreted 

1 This basic property of Regge trajectories was derived m, before the advent of DAMA. For a review of 
general properties of Regge trajectories see Ref. EU 

2 An alternative could be: a(t) ~ — ' M L (23| . 











as a barrier where the string breaks (its tension vanishes as £ —> R, producing new 
particles instead of heavier resonances, see also 123)- 

The number of thresholds is model-dependent: while the lightest one gives the domi¬ 
nant contribution to the imaginary part, those heavier promote the rise of the real part, 
terminating at the heaviest threshold. 

To illustrate the aforesaid, a toy model can be constructed from the sum of two thresholds 

a(s) = a(0) - «i(a/£i - t - y/t[) - a 2 (yjt 2 -t - y/t^), (4) 

where y/t[ is the lightest one allowed by quantum numbers, i.e. the two-pion threshold with 
t\ = 4 m 2 , corresponding to a loop diagram it the t channel |2Sj. The heavy threshold is 
chosen phenomenologically: by setting t 2 = AM 2 = lQGeV 2 , we impose an upper bound 
on the highest mass (slightly below 2 GeV) and spin (J=5) resonance lying on the given 
trajectory. The parameters ot\ = 0.6 GeV -1 and a 2 = §.5GeV~ l here were chosen such 
as to match to slope of the linear trajectory a(t) = 0.5 + t, see Fig. 1. Notice that M 




Figure 1: Left panel: typical behavior of an analytic trajectory (4) (of its real, read (r) line, 
and imaginary, green (g) line, parts) compared with a linear one, a(t) = 0.5 + t. Right panel: 
Form factor modulus (16) saturated by a single p meson pole (blue dotted line) and the p 
trajectory (4), red solid line. The form of the latter does not change when summation in 
(12) is extended to whatever large j, see Sec. 3. 


does not correspond to any physical resonance; rather it is a parameter to be fitted to 
the resonances’ spectra, as well as to the scattering data nsi- The construction of Regge 
trajectories satisfying theoretical constraints on the threshold- and asymptotic behavior, yet 
compatible with the experimental data, is a highly nontrivial problem. Simple models, like 
(4), may be helpful as a guide in a semi-quantitative analysis, as in Sec. 6. In a more 
rigorous approach of Ref. na the real and imaginary parts of the trajectories were related 
by a dispersion relation combined with the unitarity constraints on the threshold behavior 
and fits to the resonances’ masses and decay widths. 
















In the limit |u|, \t\ —> oo, u/t = const , DAMA ([TJ) scales iff its trajectories have logarithmic 
asymptotics. The simplest trajectory that combines the nearly linear behaviour at small t 
with a square-root threshold and logarithmic asymptotics is 

a(t)=a(0)- 7 ln (1±M°EP). (5) 

More thresholds (introducing more parameters, however) can be added. Asymptotically, 
a{t) ~ — ?ln(— t). The asymptotic (scaling) limit of DAMA can be easily calculated j27j 

by the saddle-point method, the saddle point, for asymptotically logarithmic trajectories, 
being located at 

= 1 / 2 . 

a(u) + a{t) 

In this limit, the ( ut ) term of the amplitude (1) goes like [27] 

D(u,t) ~ (uf)- 7ln(29)/2 . (6) 


The power in © can be fixed by the quark counting rules |29j . by which 

ln(2g) = 2 n — 1 , 

where n is the number of constituents in a collision. These numbers should not be taken 
literally (more details can be found in |2B]) since they may have more relevance to the leading, 
vacuum trajectory, while in GPD to be discussed below the main contribution comes from 
subleading trajectories. 

An interesting link between the fixed scattering scattering angle regime of DAMA 
and composite particle models in the infinite momentum frame [32J was established by 
M. Schmidt in Ref. 311 3 - 

In the model of Gunion, Brodsky and Blankenbecler (GBB) j32], the scattering amplitude 
corresponding to graph shown in Fig. 2, is given by 

A(u,t)= I dx [ d 2 k ± A'ip 1 (k ± )ip 2 (k ± + 

J o x 2 {l-x) 2 ,J 

(1 - x)q _l - xr±)ip 3 (k± + (1 - x)q ± )i/j 4 (k ± + xr ± ), (7) 

where the transverse vectors r± and q± satisfy the conditions r±-q± = 0, u = —r\ , t = — q±, 
in the infinite momentum frame 


Pi 








3 Schmidt ESI uses trajectories with a constant asymptotic limit ln(— t) —> const, relying on simple 
arguments of the wide-angle Regge behaviour in s a ( t )^ const . L ess trivial and more relevant logarithmic 
trajectories, required by the fixed angle behavior in DAMA, appeared later [77| (see also E3D- Another 
difference between the ansatze used in ESI and here is the appearance of the constant go (for more details 
on the homotopies and the role of go see m)- 











Figure 2: Diagram for (ut) channel elastic scattering, corresponding to eq. (7) 


etc., P —> +cx), and ^ are quark bound state wave functions |32j . 


The form factor F(q 2 ) and the DIS structure function Fo(x) in terms of the wave functions 
il> are ^ 

F d t f 

d 2 k± , ijj(k±)'ijj(k± + (1 - x)q±) 


F(q 2 ) = 


! 0 x(l-x) 


and 


F 2 (x) = / d 2 k ± 


2 . ^ 


1 — X 


In the limit s ~ — u = r\ and fixed t = —q\ 


A(u, t ) = 


dx 


Mrlx) / d 2 k± 


'ip(k ± )'ip(k ± + (1 -x)q_ l) 


/o x 2 (l — x) 2 ’ ,/ x(l — x) 

follows[30], where g{r\, x) ~ ijj(xr±), while for large u and fixed t one gets 2D) 

F dx 


A(u, t) 


I o z(l-x) 


girDx-^-F) f(t( l-x)). 


The parallels between dual 0 and composite models offer at least three lessons: 

1) Apart from ’’soft” collisions of hadrons (strings) [T2Jj, DAMA implicitly contains also 
the dynamics of hard scattering of the constituents; 









2 ) the variables appearing in combinations like t(l — x), should be used in constructing 
^-dependent parton distributions; 

3) Regge trajectories are non-linear, complex functions with well-defined constrains. 


3 Form factors; analyticity and unitarity 


In this section we present general properties of the form factors with emphasis on their 


analytic properties and connection with Regge trajectories, to be utilized in subsequent 
sections. 

There are various choices for the nucleon electromognetic form factors (ff), such as 
the Dirac and Pauli ff, Ff(t), Fj n (£) and F%(t), F£(t), the Sachs electric and magnetic 

ff, Fg(f), Fp{t) and Ff f (t), F£ r (t), or isoscalar and isovector electric and magnetic ff, 
FJ(f), F%(t) and F‘ M (t) , Fjtf(t), where t = —Q 2 is the squared momentum transfer of the 
virtual photon [31]. 

The Dirac and Pauli form factors are obtained from a decomposition of the matrix el¬ 
ement of the electromagnetic (e.m.) current in linearly independent covariants made of 
four-momenta, 7 matrices and Dirac bispinors as follows 


< N\J^ m '\N >= eu(p')[^F^(t) + ~P)vF?(t)]u(p), 


where m is the nucleon mass. Electric and magnetic ff, on the other hand, are suitable in 
extracting them from the experiment: 


da lab (e~N -> e~N) 


dil 

Q&m. COS 2 (fl/2)_1_ 

4F 2 sin 4 (0/2) 1 + (^f) sin 2 (0/2) 



where cc e . m . = 1/137, E is the incident electron energy, and 



or 



where p c .m. is the proton momentum in the c.m. system. 


The four independent sets of form factors are related by 



G r M (t) = G v M (t) + Gl,(t) = Ff(t) + J=f(t) = {FUt) + F?(t )] + [*?(«) + *?(«)], (9) 












Gl(t) = G‘ E (t) - Gl(t) = F?(t) + r”F“(() = |f?(t) - J?(i)] + r“[F 2 *(f) - F’(t)}, (10) 


G n M (t ) = GmW - G”„(«) = F"(t) + F 2 "(i) = [f?(f) - F”(i)] + K(t) - F*(i)), (11) 


where 


4m 2 , x 

p(n) 


. They satisfy the normalization conditions 


G p e (0) = 1; G P M ( 0) = 1 + ii p - G n R ( 0) = 0; G n M ( 0) = p n - 

^e(O) = G V e ( 0) = G^f(O) = -(1 + Pp + /in); ^7^(0) = 2 (■*■ + ~ /bOl 

Fm = 1; Tf(0) = /i p ; i?(0); F”(0) = p n] 

(o) = i?(0) = F 2 s ( 0) = i(/ip + /in); i?(0) = (/ip - /i„), 

where p p and p n are the proton and neutron anomalous magnetic moments, respectively. 

A basic ingredient of the existing models of form factors [33] is the dominance of the p 
meson pole, resulting e.g. for the isovector meson form factor to the expression 


G v (t) 


9p 


1 — t/m 


2 ’ 
P 


where g p is a constant proportional to the product of yp and 7 NN couplings. 
The next step is to include [3S other vector mesons, as well as their excita¬ 
tions, such as the isovector: p(770), p'(1450), p"(1700) and isoscalar resonances 
o;(782), u' 1420), u/'(1600), 0(1020), 0(1680) found in the Review of Particle Physics. 


The use of the trajectories implies a single ’’dual” variable instead of the parameters of 
individual resonances. The approach of Ref. PH combines the concept of Regge trajectories 
with analyticity, unitarity and resonance data analysis. An economic way to account for 
the exciting states is to use Regge trajectories, advocated in the present paper and bringing 
us close to dual models. Soon after the discovery of Veneziano’s dual model pa. attempts 
were made [3H1 to apply it to form factors. Assuming an infinity number of neutral vector 
mesons with the sequence of squared masses m 2 (n ) = m q + nm \, n — 0 , 1 ,... the following 
expression proton magnetic form factor was derived 

G Mp (t) r[l-q(t)]r[c-q( 0 )] 

p p r[i — a(o)]r[c — «(/)] ’ 

where a(t) is the p trajectory and c = 3.27. 

Unitarity constrains the threshold behavior of the Regge trajectories El (see also PH 
and references therein) as well as that of the form factors [323J 

$F w (t) ~ (t - toy-^Sait) ~(t- t 0 ) 3/2 , 

t—fto 


where 1 0 = Ait 2 . Any finite width of the resonances requires an imaginary part to be added, 
as for example was done in Ref. El 


Gm,M 

Pp 


Y a i (t — m 2 + r)i\Jb 2 — f) , 








where 7 j are resonances’ widths and the fitted values of the parameters bi are: b p = 0.28 
GeV, b u = 0.42 GeV, b$ = 0.99 GeV. This approach is close in spirit to that based on dual 
analytic model, introduced in the previous section, and illustrated in Figs. 1, showing the 
calculated modulus of the form factor resulting from a single p pole contribution (dotted 
line) and from a sequence of poles generated by trajectory (4) 


Gm p 

f- l p 


(0.5) J 

h U - a(fi)) 


( 12 ) 


This sum is similr to the pole decomposition of the dual amplitude, namely it is a sum 
of ’’reggeized” implied by duality) Breit-Wigner resonances (cf. (2)). The upper limit of 
summation includes the highest resonance lying on the p trajectory (4), however the large-|i| 
behaviour of the form factor, Fig. 3, is not affected by higher spin values (here, j > 5, up 
to infinity), from where the real part of the trajectory does not contribute any more. The 
inclusion of a large (infinite) number of poles appearing on the second sheet is important in 
dual models CHI The non-appearance of higher resonances and the transition to a smooth 
continuum can result either from an upper bound on the real part of the trajectory, as is 
DAMA [2], or from the rapid rise of its imaginary part, making the resonances unobservable. 

Similarly to the Veneziano model, form factors were derived j38] from DAMA (1): 

F(t)*= C dxx~ a ^ l ~ x) \l~x)~ 1+n , (13) 

Jo 

where n is an integer providing the correct (according to quark counting rules) large -1 1 
behavior of the form factor. 


Calculations of form factors from GPD, Sec. 4-6, are more involved and less predictable, 
especially in the case of complex Regge trajectories. 


4 Generalized parton distributions and form factors 

Form factors are related to generalized parton distributions (GPD) by the standard sum 

rules | 21 E 1 

F \(t) = / dxH q (x : £, t), 

F i(t)=J dxH q (pc, £, t). 

The integration region can be reduced to positive values of x, 0 < x < 1 by the following 
combination of non-forward parton densities dm 


H q {x, t ) = H q (x , 0, t) + H q (-x, 0, t), 
£ q (x , t) = E q (x , 0, t) + E q (—x, 0, t), 






providing 


F q (t) — / dxH q (x,t ), 
Jo 


(14) 


F q (t) = / dx£ q (x,t). 
Jo 


(15) 


The proton and neutron Dirac form factor are defined as 

F?(t) = e u F?(t)+e d F*(t ), 

F?(t)=e u F*(t)+e d F?(t), 

where e u = 2/3 and e d = —1/3 are the relevant quark electric charges. 

In the limit t — > 0 the functions H q (x,t ) reduce to usual valon quark densities in the 
proton: 

7T u ((r, t = 0 ) = u v (x), 7 i d (x, t = 0 ) = d v {x) 

with the integrals 

/ u v (x)dx = 2, / d v (x)dx = 1 

Jo Jo 

normalized to the number of u and d valence quarks in the proton. 

Contrary to H, the ’’magnetic” densities £ q (x, t — 0) = £ q (x) cannot be directly expressed 
in term of the known parton distributions, however their normalization integrals 

[ £ q {x)dx = k q 


are constrained by the requirement that the values F%(t = 0) and F%(t = 0) are equal 
to the magnetic moments of the proton and neutron, whence k u = 2 k p + k n ~ 1.673 and 
k u = k p + 2 k n ~ —2.033 follows jTT] . Explicit parameterizations for the forward and non¬ 
forward structure functions as will be presented in the next section. 


The Fourier-Bessel integral 


1 f 00 

q(x,b) = — / dv /Z t^o(^v /Z t)7-f(x,t) 
Jo 


(16) 


provides a mixed representation of longitudinal momentum and transverse position in the 
infinite-momentum frame ®E! 


5 Modelling non-forward parton distributions; 
connection with Regge-dual models 


The simplest model for the proton non-forward parton density is a factorized form 

Ti(x,t) = q v (x)Fi(t), (17) 

here q v (x) is the parton density and iq(t) is the proton form factor. It trivially reproduces 
Ff(t) and q v (x) in the forward limit, but it conflicts both with Regge (R) behavior 

n R (x,t)~ x - a ®, (18) 

valid at small x, and with the light-cone formalism, suggesting a Gaussian (G) parametriza- 
tion for non-forward parton densities mm 

n q G (x,t) = q v {x)e-^- x)t/ixX \ 

where the scale A 2 characterizes the average transverse momentum of the valence quarks 
in the nucleon. To satisfy the Drell-Yan-West (DYW) relation JflJ H2] between the x —> 1 
behavior of the structure functions and the t— dependence of the elastic form factors the 
above expression should be modified e.g. as [HJ 

H q (x,t) =q v (x)x- a '^- x) \ (19) 

where a' is the slope of a Regge trajectory (other modifications of the parton distributions 
as well as their relation to the light-cone wave function of a composite system are discussed 
in M E!)- Noticing the similarity between this expression and the relevant factor in the 
integrand of (1), we suggest the following parametrization for the t— dependent GPD: 

n q G (x,t) = q^x/go)-^ 1 -^ = (x/goy^-^fix), ( 20 ) 

where a(t) is the t— dependent part of the Regge trajectory, g 0 > 1 (x 0 in Ref. JI5]) is a 
parameter defined in Sec. 2 and f{x) is the large-x factor of the parton distribution (see 
below). For linear trajectories and go = 1, eq. (20) reduces to (19). 

Regge behavior ~ £-“(°) Q f DIS structure functions and relevant parton distributions at 
small x is well established for small and moderate virtualities Q 2 \ at higher Q 2 it is replaced 
by QCD evolution (see e.g. and references therein). The values of the ’’Regge-intercepts” 
in the parton distributions may depend on the flavor of the relevant quark. For example, in 
the global fits of MRST2002 , J|]: 

u v (x) = 0.262x _ 0,69 (1 - x) 3 - 50 (l + 3.83x°- 5 + 37.65a:), (21) 

d v (x) = 0.061x - °' 65 (1 - x) 4 ' 03 (l + 49.65x°' 5 + 8.65x), (22) 

implying «(0) = 0.69 in the w-quark distribution and 0.65 in the d-quark distribution, which 
means slightly different trajectories exchanged in the t-channel of the (Active) up and dp 
scattering amplitude, once the GPD is associated with the imaginary part of the u— channel 
quark-proton scattering amplitude. We use these expressions in our calculations below. 


The ’’magnetic” densities £ q (x,t) enter in F 2 (t) and contain new information about the 
nucleon structure, however they cannot be directly expressed in terms of any known parton 
distribution. Following m, we write them in the form £ q (x,t ) = £ q (x)x t 1 a' E = 

1.158GeV~ 2 , similar to H q (x,1) but with an extra large-x factor, i.e. 

£ u (x) = -^-(1 - x) Vu u v (x) (23) 

and 

S d (x) = -h(i ^ x)-*d v (x), (24) 

with r] u = 1.52 and g d = 0.31 fitted m to the data. The remaining constants are fixed by 
normalization: 

N u = I dx(l — xY'^uix) = 1.53; 

Jo 

N d = f dx( 1 — a:)°' 31 d(a:) = 0.82, 

Jo 

whence k u ~ 1.673, k d « —2.033. 

6 Numerical estimates 


Below we illustrate how the nonlinearity (complexity) of the trajectories and the introduction 
of g 0 > 1 affect the behavior of the calculated form factors and quark distributions as 
functions of the impact parameter b and of the Bjorken variable x. 

We first explore the effects coming from non-linear Regge trajectories by calculating 
several observable form factors and their ratios using the expressions for the t— dependent 
GPD, eq. (20). As a reference frame, we use fits from the paper [[T3j nicely reproducing 
the data. In that paper linear trajectories with the slopes oJ x = 1.098 GeV 2 in H(x,t) 
and a' 2 = 1.158 GeV 2 in £(x,t) were used. The relevant intercepts come from the parton 
distributions (21), (22), and they are equal to 0.69 and 0.65 in the u quark and d quark 
distributions, respectively M- With these parameters (and go = 0, matching the model of 
ref. jHj) one reproduces the results of Ref. J3], some of them shown in Fig.3 in thin black 
line. Next we explore the effects coming from the nonlinearity of the Regge trajectories, as 
well as from g 0 > 1, both introduced in Sec.2. By keeping the less known £ q (x,t) intact we 
vary the trajectory in 7i u (x,t) and the value of g 0 . The effect proves to be very large. For 
example, even a minor (with respect e.g. to (4)) deviation from the linear trajectory, 

a(t) = 1.026 1 - 0.02(^4m2 - t - 2 m n ), (25) 

with the forward slope a' ~ 1.098GeV~ 2 matching the linear trajectory fitted in M to the 
data, gives a sizable effect, augmenting e.g. the ratio G P M /i~i p Gd, as seen in Fig.3 (red line). 
The effect from go > 1, combined with the nonlinear trajectory, is shown in Fig. 3 (blue 
line). In general, it compensates the rise due to the nonlinearity of the trajectory, although 
the interplay of these two effects is much more complicated. The use of analytic trajectories 
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Figure 3: Left panel: Proton magnetic form factor relative to the dipole form factor (left 
upper panel) and the ratio of the magnetic to electric form factors (left lower panel). The 
black thin curves correspond to the linear trajectory of Ref. [Tlj . the red (r) curve to the 
nonlinear trajectory (24) with g 0 = 1, and the blue (b) curve to the same nonlinear trajectory 
but with g = 1.05. The data for the proton magnetic form factor G P M are from [43] (open 
squares), [H] (open circles), |4Hj (solid stars), 4 HJ (open stars), paQ) (solid circles), JH] (solid 
squares), and those for the ratio G P E /G P M are from E2J (solid circles), jJBJ (open triangles) 
and '54] (solid triangles). Right panel : Ratio of the neutron magnetic form factor to the 
dipole form factor (right upper panel) and neutron electric form factor (right lower panel) 
with curve conventions as on the left panel. The data for the neutron magnetic form factor 
G m are from [5HJ (open circles), [33] (solid circles), jHZj (open triangles), [3B] (solid triangles), 
P3J (solid squares) and [Sit (solid squares). The data for the neutron electric form factor 
are from MAMI [3T] (triangles), NIKHEF |62J (solid squares) and JLab [63; (solid circles) 
and m (°P en squares). 


and/or > 1 i n £ q (x,t) will make the situation much more complicated but, but at the 
same time, interesting. The use of a ’’realistic” trajectory like eq. (4) changes the behavior 
of the observables dramatically, requiring a complete rearrangement of the model or, at least, 
of its parameters. 

Fig. 4 shows the u quark impact parameter distribution calculated from eqs. (16), (20) 
at three fixed values of x = 0.1, 0.5 and 0.7 and for three representative trajectories: linear, 
a(t) = 1.098f (red (r) curves), ’’square root”, eq. (4) (green (g) curves) and logarithmic, eq. 
(5) with 7 = 3, (3 — 0366GeE~ 2 and t 0 = 4m 2 (blue (b) curves). The calculated distributions 
depend dramatically on the form of the trajectories, especially near the endpoints x = 0 and 
1 . 


The role of the parameter go, combined with the variation of the trajectories, can be 
seen also in Fig. 5, where the u quark distribution is plotted against x for several fixed 
values of t and of go and two typical trajectories: a linear one, a(t) = 1.098f, and Q. Here 
the dependence on the form of the trajectories is less pronounced then in Fig. 4, where 
integration in t is involved. 






















Figure 4: Impact-parameter distribution of the u quark matter at x — 0.1 (left panel), 
x = 0.5 (middle) and x = 0.7 (right panel). The curves on the left panel correspond (from 
the top) to calculations with: a(t) = 1.098f (red, solid line); eq. (4) (green, dotted line) and 
the logarithmic trajectory eq. (5) with 7 = 0.3, f3 = 0.366 and to = 4m 2 (blue, dotted line). 
Curves on the middle and right panel correspond to calculations with eq. (4) (green dotted 
line) and the above linear trajectory (red solid line). Those with the logarithmic trajectory 
(5) for x > 0.3 go off the common trend. 

u u 




Figure 5: Longitudinal momentum distribution of the u quark matter calculated with a 
linear trajectory a(t) = 1.098f (left panel) and eq. (4) (right panel). The curves (from the 
top) correspond to: [t = —0.7GeV 2 , go = 1] (red line), [t = — 0.9GeV 2 , go = 2 ] (blue line), 
[t = — 0.7GeV 2 , g 0 = 5] (green line), [t = —3 GeV 2 ,g 0 = 1] (pink line). 


Given the available freedom in the choice of the trajectories, the present calculations can 
serve only as an indication of the existing trends, rather then fits or predictions to the data. 


7 Conclusions and outlook 


With this paper we wish to emphasize the important role of the analytic properties of the 
strong interaction theory, manifest here in the form of the Regge trajectories. The wide¬ 
spread prejudice that the trajectories are linear has different sources: 1 ) the masses of the 
resonances he on approximately linear trajectories; 2) the Veneziano and string models [H9] 
provide a theoretical basis in favor of this behavior; 3) relevant calculations are simple. On 
the other hand, the theory demands that the trajectories be analytic functions of their argu- 



















merits with threshold singularities imposed by unitarity, and that the asymptotic behavior be 
compatible with the polynomial boundedness of the amplitude. Trajectories satisfying these 
constraints and fitting the data both for positive (particles spectra) and negative (scattering 
data) values of their arguments are known from the literature (see, e.g., [ITB1 fTT|i. In this 
paper we show that they affect considerably the calculated GPD and their moments. More 
work is needed to specify and quantify the role of various (parent and daughter) subleading 
trajectories (of poles and cuts). 

Another message of this paper is that, in a certain kinematical region, the integrand 
of the dual amplitude with Mandelstam analyticity (DAMA) (1) can be identified with a 
generalized parton distribution (GPD), the integration variable being associated with the 
parton longitudinal momentum, as suggested in For fixed Q 2 and s, the integrand of 
DAMA (1) (a GPD?) has the form 

(i/9o)- aW1 '-»-”((l - n/9o) m , (26) 

where go > 1, and n and m are reaction-dependent constants. The first factor in (26) is the 
sniall-x Regge-behaved term, while the second one is the familiar large-x term. An immedi¬ 
ate observation is the similarity between the moments of the GPD (form factors), eq. (14) 
and expression (13), derived from DAMA. I wonder if this could mean a bootstrap relation 
between the antiquark-hadron scattering amplitude (GPD) and the hadron-hadron ampli¬ 
tude resulting from integration (1) of a GPD ? The appearance of x in the exponents of the 
(modified) structure functions, (or parton distributions) may have interesting consequences 
by itself - both for theory and phenomenology. In any case, a better understanding of the 
physical meaning of the variables appearing in GPD, especially with non-zero skewedness, 
£ 7 ^ 0 , is needed. 

The approach in this paper combines elements of the analytic S— matrix theory, namely 
Regge poles and duality, known to be efficient for ’’soft” collisions at large distances, with 
the small-distance partonic picture. The interface and merge of these seemingly orthogonal 
approaches may bring new ideas about the transition form perturbative to non-perturbative 
physics. Much of this information is encoded in the form of the complex Regge trajecto¬ 
ries. Recently, explicit models for deeply virtual Compton scattering amplitudes (DVCS) 
appeared in the literature (Ii|. Their imaginary part can provide additional information 
about GPD. 

In a perspective one can think of extending the asymptotic Regge pole model to the 
low-energy resonance region by incorporating a dual amplitude, e.g., DAMA. The crossing- 
symmetric properties of dual amplitudes will make possible the inclusion of the t > 0 region 
in DVCS and resulting GPD. It may also help to connect the high-energy (Regge) region 
with the low-energy (resonance) domain, where new data from JLab is expected. 
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